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We present a study of the dynamic structure factor of the antiferromagnetic spin- 1/2 Heisenberg 
chain at finite temperatures and finite magnetic fields. Using Quantum-Monte-Carlo based on the 
stochastic series expansion and Maximum-Entropy methods we evaluate the longitudinal and the 
transverse dynamic structure factor from vanishing magnetic fields up to and above the threshold 
B c for ferromagnetic saturation, as well as for high and for intermediate temperatures. We study the 
field-induced redistribution of spectral weight contrasting longitudinal versus transverse excitations. 
At finite fields below saturation incommensurate low-energy modes are found consistent with zero 
temperature Bethe-Ansatz. The crossover between the field induced ferromagnet above B c and the 
Luttinger liquid below B c is analyzed in terms of the transverse spin-dynamics. Evaluating sum- 
rules we assess the quality of the analytic continuation and demonstrate excellent consistency of the 
Maximum-Entropy results. 
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I. INTRODUCTION 

The antiferromagnetic spin- 1/2 Heisenberg chain 
(AFHC) is one of the most intensively studied strongly 
correlated quantum many body systems. In the pres- 
ence of an external magnetic field, its generalization to 
anisotropic exchange, the XXZ model, reads 

H = J$>f Sf +1 + j(S+S w + S^S+ +1 ) BSf], (1) 
i 

where J is the antiferromagnetic exchange interaction 
with an anisotropy ratio A, S^' z arc the spin opera- 
tors on site I of the chain, and B = g[i B hH is the 
magnetic field. From a materials perspective SrCuO^ 
(J/k B « 260UFO, Sr2Cu0^iM (J/k B « 22001T), and 
Cu(C 4 H 4 N2)(N03)2^ (J/k B « 10.7K) are topical exam- 
ples of both, low- and high- J AFHC compounds which 
have been studied intensively. Recently, dynamical cor- 
relation functions of the AFHC have become accessible to 
a variety of high resolution spectroscopics at finite tem- 
perature and in the presence of external magnetic fields, 
eg. inelastic neutron scattering (INS) 7 -, high-field nuclear 
magnetic resonance (NMR)^ ' 10 ' 11 , muon spin-resonance 
(|tiSR)i2, and magnetic transpor t 13 ' 14 . 

The AFHC is integrable, including the cases of h ^ 
and A / 1. Bethe-Ansatz (BA)i^ has been used to 
investigated its ground state properties. Static ther- 
modynamic quantities, eg. the specific heat, the mag- 
netic susceptibility and the magnetization have been 
investigated by several methods including thermody- 
namic BA, Quantum Monte Carlo (QMC), as well 
as transfer-matrix and density-matrix rcnormalization 
(DMRG) group, se o 17 ' 18 ' 19 and refs. therein. 

Fractionalization of the spin excitations into multi 
spinon states is a fingerprint of the AFH C 20 : 21 . At zero 



temperature, T = 0, numerical analysis of these excita- 
tions has been carried out in many studies using exact di- 
agonalization (ED) of finite AFHCs, see eg. refs. 22 : 23 ' 24 , 
including the effects of h and A, as well as by dynami- 
cal variants the DMR G 25 ' 26 . In principle, also BA allows 
to determine dynamical correlation functions, however 
calculating the corresponding matrix elements is highly 
non-trivial and progress has been made only recently. By 
now analytic expressions for dynamical spin correlation 
functions are available for the two 27 i 28 i 29 and the four- 
spinon sector^ii 3 ^ 3 . at A = 1, h = and T = 0. In ad- 
dition, determinant approache a 34 ' 35 allow for numerical 
treatment of two 36 i 37 : 38 and many-spino n 39 : 40 states of 
the XXZ chain in finite magnetic fields, at T = 0. Finally, 
mapping to field theory in the continuum limilj 41 - has been 
used to study the small-g behavior of longitudinal dy- 
namical structure factor in the gapless regim e 42 : 43 ' 44 . 

At finite temperatures, the dynamical correlations 
functions of the AFHC remain an open issue. The dy- 
namical structure factor S a/3 (q,ui) has been studied by 
complete ED of small system a 45 ' 46 in the context of spin 
diffusion, se o 13 ' 47 . and refs. therein. However, such 
analysis is limited by finite size effects to k B T > J. 
Recently, finite temperature real-time auto- and next- 
nearest neighbor correlation functions have been accessed 
by DMRG method a 48 ' 49 . However, the time range of such 
calculations is limited, as the spectrum of the reduced 
density matrix used to truncate the Hilbcrt space be- 
comes dense. In this respect QMC remains a key tool to 
evaluate the S a/3 (q,cu), for system sizes which are close 
to the thermodynamic limit, over the complete Brillouin 
zone, and at finite temperatures, with the limitations set 
primarily by the analytic continuation of imaginary-time 
dataS -. QMC analysis of S af3 (q, u) has been carried out 
for h = O 51 ' 52 ' 53 , results for h ^ 0, however, are lacking. 

The purpose of this work is to shed more light 



2 



the finite-temperature dynamical structure factor of the 
AFHC in the presence of external magnetic fields using 
QMC. The paper is organized as follows. In section ITT1 we 
briefly summarize the QMC approach we use. Section 
IIIII we analyze the transverse and longitudinal structure 
factor versus temperature and magnetic field. In section 
IIVI we consider several sum rules of the AHFC. We sum- 
marize and conclude in section fVl 

II. METHODS 

In this paper we present results for the Fourier trans- 
form of the dynamic structure factor 

/oo 
dte iut S a0 (q,t) , (2) 
-OO 

where t refers to real time and S a/3 (q,t) = (S%(t)St q ) 
with S q = e~ lqr being the spin component a at 
momentum q. In the following we discuss both, the longi- 
tudinal and the transverse structure factor, i.e. a/3 = zz 
and a/3 = xx. These components are directly accessible 
to unpolarized neutron scattering. Other types of trans- 
verse components, eg. a/3 — H — require polarized neu- 
trons and will not be considered here^. For a/3 = xx(zz) 
S a @(q,us) is related to the imaginary time structure fac- 
tor S al3 (q,T) = (Sg(T)S_ q ) through the integral trans- 
form 

S af3 (q,T) = - / du>K(u>,T)S af3 (q,cj) (3) 
Jo 

with the kernel K{uj,t) = e -Tt " + e ~^~ T ^ . 

The imaginary time correlation functions S a ^(q, r) will 
be evaluated by QMC. Here we employ the stochastic se- 
ries expansion (SSE) introduced by Sandvik et al£H. This 
method is based on a particular form of the high tem- 
perature series expansion of the partition function which 
can be sampled efficiently 56 . To this end, the Heiscnbcrg 
model is rewritten in terms of N b bond operators 

N„ 

H=-jJ2(Hi,b + H 2 , b ) (4) 

6=0 

where H hb = 1/2 - S Z S^ +1 and H 2 , b = (SfS^ + 
S^f S^ +1 )/2 refer to the diagonal and off-diagonal parts 
of H within the S z basis. Within this notation the par- 
tition function is 

*=£££^h]W>> (5) 

a n S n k=l 

where \a) = \Sf , . . . , S^) refers to the S z basis and S n 

S n = [a 1 ,b 1 ][a2,b2]...[a n ,b n ] (6) 

is an index for the operator string Y\Z = iH akbk , labeling 
each specific product of operators where £ {1, 2} and 
b k e {!,... N b }. 



The operator string is subject to importance sampling 
by a Metropolis scheme, splitted into two different types 
of updates: a diagonal update which changes the number 
of diagonal operators H\ bk in the operator string and a 
cluster type update which performs changes of the type 
H\.b k «-> H2,b k - On bipartite lattices the latter, so-called 
loop update guarantees an even number of off-diagonal 
operators H2.b k in the expansion. This ensures positivc- 
ness of the transition probabilities. 

Imaginary time spin-correlation functions S"j(r) = 

(e rH S«e- TH sP), for r G [0,0), and lattice sites i,j can 
also be sampled by the SSE using that^l 

\m=Q ' v ' I w 

where, in contrast to the thermal average (...}, the brack- 
ets (. . .)w with subscript W refer to averaging over oper- 
ator string configurations with weights generated by the 

Metropolis scheme. The quantity (7,- • (to) in eqn. (J7J) is 
the static real space correlation function 

M = — T £ 5f(p)5f (p + to) (8) 

' J 71+1 ' — ' J 

p=0 

which, in the case a/3 = zz, can be measured within 
the diagonal or slightly more efficient within the loop 
update^ where m,p < n refer to positions within the 
operator string and Sf (p) refers at the intermediate state 
\a(p)) = rife=i Ha k ,b k \ a ) of the expansion. Transverse 
correlations a{3 = xx, H — , — h in a code working in the 
S z basis can only be accessed by following the loop within 
the extended configuration space (see§£). Finally, the 
Fourier transform S al3 {q, r) = £\ e~ lqj S^{t) / N is used 
as input for the inversion problem eqn. (|3"|). 

Extracting S a " (q, u) from eqn. ([3]) is notoriously com- 
plicated because of the QMC noise and the singular na- 
ture of K(u>,t). However, this problem can be handled 
successfully by Maximum Entropy methods. Here we 
use Bryan's algorithm^ which is specifically designed for 
over-sampled data sets and therefore well-suited to treat 
QMC results. For details we refer the reader to the Ap- 
pendix. 

III. RESULTS 

In this chapter we will present results for the trans- 
verse and longitudinal dynamic structure factor at finite 
temperatures, in the range of T — J/20 . . . J and mag- 
netic fields below and above the saturation field B c . All 
QMC calculations refer to systems with 128 sites, typi- 
cally with one billion Monte-Carlo updates (one diagonal 
and sufficient^ loop updates), distributed over 1000 bins. 
Only 50-100 T-points were extracted for each tempera- 
ture in order to prevent over-sampling of the relatively 
short expansion orders at elevated temperatures close to 
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FIG. 1: 3D plot of the longitudinal dynamic structure factor by QMC + MaxEnt as function of frequency to and wave vector 
q. Temperatures and magnetic field in units J: a) T=l and B=0, b) T=0.25 and B=0, c) T=l and B=l, d) T=0.25 and B=l, 
e) T=l and B=2 and f) T=0.25 and B=2. 



T = J. An indication for over-sampling is given by diag- diffusion is a long-standing issue with no final answer as 
onalizing the covariance matrix which exhibits vanishing of today. For a recent review on the present status and 
eigenvalues in case of statistically dependent data. related referencs we refer to^. 

Unfortunately QMC is too sensitive to the default 
model for the MaxEnt continuation in the small-g, w 
A. Longitudinal Dynamic Structure Factor S zz (q,u)) regime^ to elucidate the issue of spin-diffusion. Yet, we 

would like to mention agreement of our results regard- 
ing the frequency-transformed autocorrelation function 
Sq z (uj), i.e. the q-integratcd dynamic structure factor 
(not shown within this work) with previous QMC, per- 
formed at B = 0, high temperature series expansion^ 
and TMRGi2. These results exhibit a w -o.3...-o.4_ 
divergent behaviour which bears resemblance to the 
phenomcnological approaches by Blocmbcrgcn^ and de 
Gennes^ who predicted lj" 1 / 2 . 

Next we consider lower temperatures, i.e. T = J/4. 
As is obvious from Figs. QJ>) and|2fr), spectral weight is 
removed from the long wave-length regime in this case. 
Both figures demonstrate that most of the spectral weight 
is confined within the two-spinon boundaries with how- 
ever still an appreciable intensity below the lower bound- 
ary. This is consistent with findings reported inaa. In 
contrast to T = J we find a strongly enhanced spectral 
weight at q = n owing to the increase of the antifer- 
romagnetic correlation length-^ which is consistent with 
the autocorrelation function reported in ref4^. 

In the limit (q, to) — > (it, 0) we find indications for di- 
verging behaviour of S zz (tt,lu) with decreasing tempera- 
tures. This is shown in more detail in Fig. [31 scanning 
a wide range of temperatures from T = J to T = J/20. 
As can be seen, the spectrum consists of an upturn for 
uj — > and a peak at finite uj. The latter peak shifts to 
lower energies while gaining sharpness as T — > 0. For 



In Fig. [T] and [2] we show the longitudinal dynamic 
structure factor both, as a 3D and a contour plot for 
two different temperatures T = {J, J/4} and three dif- 
ferent magnetic fields B = {0, B c /2, B c }. The solid lines 
displayed in the contour plots for B < B c refer to the 
upper and lower boundaries of the two-spinon spectrum 
as obtained from BA selection rules^. For zero mag- 
netic field they enclose a region which, within 2-spinon 
calculations, contains about 73% of the zero tempera- 
ture spectral weight—. We will now focus on each of the 
magnetic fields separately. 



1. The case B = 

At zero magnetic field and high temperatures, i.e. Figs. 
QJ.) and (2k), we find a strong broadening of spectral 
features. While the region of finite spectral weight re- 
mains bounded from above by J-tt sin |g/2|SS, significant 
weight appears below the lower two-spinon boundary 
^ sin |g| set by de Cloizeaux-Pearsor^i. Most notewor- 
thy, high spectral weight occurs for q, uj — > 0. This inten- 
sity is related to spin conservation which dominates the 
long wave-length dynamics in the quasi-classical regime 
fcgT ^> J. The question whether the long wave-length 
spin dynamics in the AFHC can be described by spin 
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FIG. 2: Contour plot of the longitudinal dynamic structure 
factor as function of frequency u and wave vector q. Tem- 
peratures and magnetic fields in units of J: a) T=l and B=0, 
b) T=0.25 and B=0, c) T=l and B=l, d) T=0.25 and B=l, 
e) T=l and B=2 and f) T=0.25 and B=2. For B < B c the 
solid lines are zero temperature excitation boundaries by the 
Miiller-ansata^ while at critical fields the exact zero temper- 
ature 1 — cos(g) dispersion^ is shown. 



T — > 0, Fig. [3] suggest that the peak will merge with the 
zero-w upturn to form a single divergence at to — > 0, as 
predicted by 2-spinon calculations at T = which lead 
to S zz (ir,uj) ~ cj _1 — . A similar peak at finite u> was 
observed also in ref.— . However, smaller systems sizes in 
that case, i.e. N = 32, render the zero-w upturn into a 
shoulder only. Biasing the default model by several sum- 
rules, it was shown in refi^, that S zz (ir, u>) on 32-site sys- 
tems could be obtained with only a single peak at finite 
oj. Recent SSE-QMC on 128-site systems at B = show 
only a single rounded maximum, centered at oj = 0^. 
While all these findings are consistent with the forma- 
tion of a zero-w divergence as T — > 0, they show that the 
details of the low-w spectrum are subject to details of 
the MaxEnt approach. Nevertheless, we will detail later 
that our results are consistent with several sum-rules, in- 
cluding those which arc particularly sensitive to the low 
frequency behavior of the spectrum. 
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FIG. 3: Zero-field dynamic structure factor at q — n for four 
different temperatures (in units of J) {f, 0.25, 0.1, 0.05}. As 
the temperature decreases we find an increased divergent be- 
haviour for oj — > as predicted by 2-spinon calculations. In 
addition there is a low frequency peak which shifts to lower 
energies while steadily gaining sharpness. Note that the data- 
set for T = 1 and for T = 0.25 was multiplied by a factor of 
eight, respectively five for illustrative reasons. 



2. The case B = B c /2 

Figs. 11121 c) and d) depict the longitudinal structure 
factor at half of the critical field. The impact of a fi- 
nite magnetic field is fourfold. First, at zero momen- 
tum the longitudinal structure factor is proportional to 
the square of the field-induced magnetization at zero fre- 
quency, i.e. S zz (q = 0, w) ~ (S z ) 2 6(oj). To focus on the 
remaining spectrum, we have chosen to skip the single 
wave vector q = in all 3D, as well as contour plots 
of S zz (q,uj) for B ^ 0. Second, longitudinal excitations 
with will have decreasing matrix elements with 

increasing magnetic field. This is consistent with the 
evolution of the overall scale in Fig. Q] a)-e) and b)-f). 
Third, longitudinal spin-excitations at the zone boundary 
are energetically unfavorable in a magnetic field. In fact, 
at low temperatures a gap can be observed at q = n, 
which is proportional to the magnetic field^ (see Fig. 
[5]d)). Finally, a soft mode occurs at an incommensu- 
rable wave vector q s = 7t(l - 2{S Z }) (see Fig. [2] d)). 
This can be understood in terms of the Jordan- Wigncr 
fermionic description of the AFH C 65 ' 66 i 67 , where S z is re- 
lated to the fermion density and the magnetic field plays 
the role of a chemical potential driving incommensura- 
bility. This finding is consistent with refi^, with inter- 
acting spin-wave calculation^ as well with finite system 
diagonalization££. The role of temperature is evident. At 
high temperatures, i.e. T = J in Fig. H]c), S zz (q 7 uj) is 
rather featureless and extends clearly beyond the bound- 
aries set by the two-spinon continuum. This changes as 
the temperature is lowered to T = J/4, Fig. H]d), where 
the spectrum is far more confined to within the dCP 
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FIG. 4: Transverse dynamic structure factor by QMC + MaxEnt as function of frequency u) and wave vector q. Temperature 
and magnetic field in units of J: a) T=l and B=l, b) T=0.25 and B=l, c) T=l and B=2, d) T=0.25 and B=2, e) T=l and 
B=2.5 and f) T=0.25 and B=2.5. 



boundaries and displays more pronounced features. In 
particular, the weight is enhanced as (<?, oj) — ► (jr, J). 



3. The case B = B c 

For B > B c and T = 0, the statistical operator of 
the AFHC is pure and corresponds to the fully polarized 
state, i.e. S zz (q,u) = N(l/4) 5 q fi5(uj). Additional finite 
spectral weight for q,oj ^ will occur only for T > 0. 
To observe this we have again removed the wave vector 
q = from Figs. [TJ H]e), f), which are at the critical 
field. Indeed, on lowering the temperature from panel e) 
to panel f) in Fig. [TJ the remaining total spectral weight 
decreases. Apart from this the higher temperature spec- 
trum is rather featureless, while the lower temperature 
spectrum clearly resembles the exact zero temperature 
dispersion of (1 — cos(q))^ (see Fig. E]f)). This exci- 
tation has a constant spectral weight 2ir/N for q ^ 0, 
which vanishes in the thermodynamic limit. 



B. Transverse dynamic structure factor S xx (q,u)) 

In Fig. [4] and [5] we show the transverse dynamic struc- 
ture factor as 3D and as contour plots for identical tem- 
peratures T/J = {1,0.25} as for the longitudinal dy- 
namic structure factor, however for a different range of 
magnetic fields B/J = {1,2,2.5}. For vanishing mag- 
netic field we refer to Figs. [TJ [2] a), b) for S xx (q,ui) 
which is identical to S zz (q,uj) at B = due to SU(2) 
invariance. 



1. The case B — B c /2 

First, we note that the results for S xx (q, oj) in Figs. 141 51 
a) and b) are clearly different from those for S zz (q,oj) in 
Figs. 11121 c) and d) at identical magnetic fields. This is 
to be expected, since the application of a finite magnetic 
field breaks the SU(2) invariance of the AFHC. Second, 
long-wave length transverse spin-excitations will experi- 
ence the Zeeman energy due to the magnetic field, which 
leads to a spin gap of size B/J at q = 0. This has to be 
contrasted against the gap at q = tt in S zz (q,oj) at finite 
fields in Figs. [U [21 d). Third, and as for the longitudi- 
nal case a field driven zero mode at q = q s can be seen 
in [2] a), b) - with a rather low intensity as T — > 0. In 
contrast to the longitudinal case, this mode develops out 
of the zone zenter and moves to the zone boundary with 
q s = 27r(5 z ^S. 

Even though it can be misleading to compare MaxEnt 
data based on different QMC data sets quantitatively 
due to the underlying different statistic quality, we no- 
tice enhanced spectral weight near the zone boundary in 
Fig. [4]b) compared to zero magnetic field in Fig. Q]b), 
which means that a weak uniform field strengthens the 
antiferromagnetic order in the transverse strucure factor. 
This effect was also observed vd& for small fields and by 
Karbach et al. within the static structure factors^. 



2. The case B > B c 

At intermediate fields selection rules^ allow for a fairly 
complex distribution of spectral weight as is also obvious 
from the solid lines in Figs. [5] a) and b). In contrast to 



6 



3 



a) 0.7 1 


(3.0 d) 


b) 1^31 


4.3 g) 


C) 0.9 


4.9 f) 



7T./2 n 



turn spin-systems recently^^^^^^^^^ 3 *^^. 

These systems feature a gapful zero-field state with the 
lowest triplet branch 'condensing' as the field is increased. 
For the AFHC, this scenario is reversed, i.e. decreas- 
ing the field through the critical value for complete po- 
larization B c , the magnons condense at q = it and the 
system switches from a gapful state to a Luttingcr liq- 
uid of deconfined spinons. Obviously, the latter docs 
not represent a true gauge-symmetry broken state, since 
(i) ID-correlation functions decay algebraically and (ii) 
the magnons above B c are constrainted by a hard-core 
repulsior&S&S 9 .. 

As the temperature is lowered, the thermal smearing 
of the approximately quadratic dispersion at q = tt for 
B = B c is reduced, see Figs. [5] c) to d). For the 
momentum-integrated structure factor this will lead to 
a critical increase of the density of states at uj = as 
temperature is lowered. Clear indications of this critical 
enhancement as T — ► have been observed in recent mea- 
surements of the transverse relaxation rate T^ 1 in nu- 
clear magnetic relaxation experiments at B = B c on cop- 
per pyrazinc dinitrate Cu(C 4 H 4 N2)(N03)2 (CuPzN)±i. 



FIG. 5: Transverse dynamic structure factor by QMC + Max- 
Ent as function of frequency u) and wave vector q in a contour 
plot. For parameter details see text or Fig. [4] The solid lines 
for half critical field a) and b) are zero temperature excitation 
boundaries of different B A selection rules (see^) . For B > B c 
the one-magnon cosine dispersion is shown. 



this, above the saturation field and at low temperatures, 
a straightforward picture emerges (see Fig. 14151 d) and 
f)). In this regime and for T — > the systems is fully 
polarized. In that case the elementary excitations are 
non-interacting ferromagnetic one-magnon states, lead- 
ing to a dispersion E(k) — Jcos(fc) + B in the trans- 
verse structure factor, with a momentum-independent 
spectral weight^. For finite T, we find that this pic- 
ture is modified in two ways. First, significant thermal 
broadening occurs, which as e.g. in Fig. [5]e), at B = 2.5 
and T = J can lead to a complete closure of the zone 
boundary spin-gap. Second, and as can be seen in Figs. 
@] d) and f), there is a substantial wave- vector depen- 
dence of the spectral weight in the cosine-signature of 
the one-magnon state. The latter is due to the elemen- 
tary one-magnon states being excitated in a polarized 
background which contains thermal fluctuations^. Fi- 
nally, we emphasize the difference in the evolution of the 
overall spectral weight, contrasting longitudinal versus 
transverse excitations. While in Fig. [T]the weight of the 
excitations decreases with increasing field, this is not so 
in Fig. |U 

Figs. 2] and [5] bear a close resemblance to the concept 
of field-induced Bose-Einstein condensation of triplets, 
which has been under intense scrutiny for several quan- 



IV. SUM RULES 

Sum rules have been used extensively for the AFHC to 
evaluate the contribution of 2- and 4-spinon excitations 
to the spectral weight of the dynamical structure factor 
at T=0 28 i 33 i 90 i 91 . For the present work sum rules can be 
applied to assess the quality of the analytic continuation. 
We will focus on the sum rules for the static structure 
factor S a P(q) and the static susceptibility x a ^{l) which 
obtained by integral transforming the dynamical struc- 
ture factor— 

1 f°° 

S a0 (q) = ~ duj(l + e~ luJt )S a P{q,Lo) (9) 
Jo 

X aP (q) = - / duu-\l - e~^)S^{q,w). (10) 

While S a,3 (q 1 Lo) on the right-hand side of cqns. © 
and (flC)]) involve MaxEnt-data, the static structure fac- 
tor S a @(q) in cqn. (J9j) is calculated from a real-space 
Fourier transformation of the equal-time correlation func- 
tions and the static susceptibility x a/3 (<l) m ecm - (PH can 
be evaluated from the Kubo integral 

X af3 (q)=J2 eiqr J o P dr(s?(r)SS(0)). (11) 

of the imaginary time QMC-data. I.e. both, S a "(q) and 
X a ^(q) are obtained from QMC-data which is indepen- 
dent from the MaxEnt continuation. In particular the 
static susceptibility should provide for a clear consistency 
check regarding the low energy features in the zero field 
dynamic structure factor at k = 7r as shown in Fig. [3] 
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FIG. 6: Comparison of transverse and longitudinal static sus- 
ceptibility/structure factor (symbols) and sum rules (lines) for 
T=J/4 and four different magnetic fields B/J = {0, 1, 2, 2.5} 
(from top to bottom). All sum rule results are within the 
error bars of the static quantities which are within symbol 
size. 



In Fig. [5] we compare the left- and right-hand sides of 



cqns. ([9]) and ([TO]) both, for the longitudinal and trans- 
verse components, i.e. a/3 = zz and a(3 = xx. First, we 
emphasize that the numerical values for S zz (0), S zz (tt), 
X 2Z (0) and x zz ( 7T ) which we have obtained at zero mag- 
netic field are consistent with those reported in refsi 53 i 64 
and corroborate the parameters of scaling relations^ 

S zz {-k) = D s \n{T s /T)§, D s = 0.094(1), T s = 18.3(5) 



X zz (n) = ^HT x /T)i, D x = 0.32(1), T x = 5.9(2) 

for T = J/4. Second, Fig. [6] proves an excellent agree- 
ment of QMC data involving analytic continuation to 
that free of the MaxEnt procedure. We have found this 
agreement for all temperatures and all fields investigated, 
including those not depicted here. All differences lie 
within the error bars of the static quantities which is 
remarkable, given that the typical MaxEnt error is esti- 
mated to be ~10-20%S. We note that we have performed 
this sum-rule check for various MaxEnt procedure, i.e. 
historic, classic, and bryan (see Appendix |Aj and found 
the same level of agreement. 



V. CONCLUSION 

In conclusion, using MaxEnt continuation of QMC re- 
sults, we have analyzed the evolution of transverse and 
longitudinal spin excitations of a AFHC with 128 sites 
at finite temperatures and magnetic fields up to and 
above the saturation field. Our results are consistent 
with and complement similar studies using small system 
ED and zero-temperature BA. In particular we have de- 
tailed the difference between longitutinal and transverse 
excitation as a function of the magnetic field and tem- 
perature. Moreover we have considered the field induced 
magnon 'condensation' at the saturation field and the 
occurence incommensurate zero-modes. These investi- 
gation may be of relevance to high-field NMR data on 
AFHC materials^ as well as to inelastic neutron scatter- 
ing experiments. 

Several open questions remain. While the issue of spin- 
diffusion has been out of reach in this work, future analy- 
sis should improve the resolution of the MaxEnt, in order 
to access the line-shapes at small q. This also pertains to 
the form of the low-w spectrum of the zero-field dynamic 
structure factor at q = tt/2. Finally it will be interesting 
to perform similar calculations for various generalizations 
of the AFHC including anisotropy and disorder. 
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APPENDIX A: MAXIMUM ENTROPY METHOD 

In this appendix, and for completeness sake, we give 
a brief account of our MaxEnt approach. In MaxEnt we 
minimize the functional 94 



1 



(Al) 



For perfectly uncorrelated QMC data % 2 ls the least- 
square difference between the data S al3 (T), with standard 
deviation o~ T , and the transform of the trial spectrum 
A(tu) to imaginary times using the Kernel K{lo, t) 



X = 



E 



1 2 



(A2) 



where for brevity we disregard the q-dependency. In prin- 
ciple imaginary-time output from the QMC is correlated 
and needs to be transformed into an eigenbasis of the 
covariance matrix prior to using eqn. (|A2jl in order to 
work with decorrelated data. However, we have observed 
that diagonalizing the covariance matrix has neglcgiblc 
impact on the spectra which we have analyzed. There- 
fore we have decided to ignore off-diagonal elements of 
the covariance matrix. 

The second term on the right-hand side of eqn. (|A2[) 
contains the Shannon entropy 



5 = E 



A(lu) - m(uj) - A(u) In 



AH 

m(uj) 



(A3) 



with respect to a default model m(uj) which prevents 
overfitting of the data. We have used a simple flat default 



model for all calculations which was itcrativcly adjusted 
to match the Oth moment of the trial spectrum. This is 
different from ref.— , where several sum-rules have been 
used to apply additional bias to A(oj). 

The choice of the Lagrange parameter a has been dis- 
cussed extensively^. So-called classical and historic ap- 
proaches use Bayesian logic to fix one a for the most 
probable spectrum A a (u>). More generally however, a 
probability distribution P \a\S a P(T)\ exists^, which de- 
termines the most likely spectrum through the average 



S a(3 (Lu) = J daA a (u)P [a\S aP (T)] . (A4) 



We have analyzed our results in term of all three ways 
to choose a. We found classic and averaged spectra to 
be identical, indicating that P [a\S a/3 (t)] is very sharp 
in our case, which supports the statistical quality of the 
underlying QMC data. As to be expected, for the historic 
approach we found somewhat smoother results with a 
tendency to under-fit the data, thus all shown results are 
based on averaged spectra. 

The minimization of (|A1|) is done via multi- 
dimensional Newtown iterations. However, following 
Bryan's work^ we have reduced the effective search di- 
rections by a singular value decomposition of the kernel 
K = UY>V T down to typically 10-20 of the largest eigen- 
values of E, depending on the temperature. This leads 
to a significant speed up of the algorithm. 
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